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AMPLIFICATION BY WAVE DISTORTION OF THE DYNAMIC 
RESPONSE OF VAPORIZATION LIMITED COMBUSTION 
by Marcus F. Heidmann 
Lewis Research Center 

SUMMARY 

The response of a combustion process to acoustic oscillations which have distorted 
wave shapes is analyzed. A combustion process where burning rate is proportional to 
an exponential power of a Reynolds number is assumed, and specific application is made 
to drop vaporization in rocket combustors. 

Analytical solutions are derived for the in-phase and out -of -phase response factors 
(components of the burning rate in-phase and out -of -phase with pressure). The response 
is shown to depend on the harmonic coefficients and phase angles in the Fourier series 
used to describe distorted wave shapes. The analytical solutions are compared with 
exact numerical evaluations. 

The addition of harmonic content to acoustic oscillations can increase the in-phase 
response factor by an order of magnitude above that for sinusoidal wave shapes. The 
largest increase occurs when the harmonic components and acoustic particle velocities 
are all in-phase. The out -of -phase response is zero for this condition, but with phase 
shifts it has a finite value. 

The effects of wave distortion on the response factors depend on the relative axial 
velocity of vaporizing drops. Large effects are observed when relative velocities are 
less than about 30 meters per second or 500 feet per second. 


INTRODUCTION 

Nonlinear instability limits of rocket engine combustors are one of the least under- 
stood properties of unstable combustion. These limits are usually established by ex- 
ploding bombs of various sizes within the combustor. The engine remains stable to a 
bomb disturbance below some threshold level, and self -sustained instability is triggered 
above that level. Although simple in concept, such stability rating has uncertainties. 



Bomb size alone is not the only variable. The type of bomb and its location also affect the 
results. In addition, instability is occasionally triggered by a bomb which normally does 
not cause instability, or it can be self -triggered by the combustion without bombs. The 
cause of instability and the variables affecting the stability limits are of vital concern 
in the design of high performance rocket engines. 

The purpose of this study is to demonstrate the importance of wave distortion to 
both the analysis of instability and to the design of stable combustors. Distorted waves 
(pressure wave shapes which depart from pure sinusoidal forms) are generated by the 
bombs that trigger instability. Also, distorted wave shapes are usually observed during 
acoustic mode instability. The distortion usually increases with the amplitude of the 
oscillations. In the limit a shock wave or other highly nonsinusoidal wave form de- 
velops. Wave distortion, however, is not limited to high amplitudes. Instability in some 
engines exhibits simultaneous resonance of several acoustic modes - usually at harmonic 
frequencies - and distorted wave shapes are seen at very low amplitudes . 

The effect of wave distortion on the acoustic amplification that can be expected from 
a drop vaporization process is examined. Steady combustion in liquid rockets is usually 
assumed to be vaporization limited, and the quantitative agreement between theory and 
experiment has been good (ref. 1). Vaporization -limited combustion is also expected for 
nonsteady combustion. Although this has been qualitatively confirmed, there are quan- 
titative deficiencies in theoretical analyses of the problem. Dynamic analyses of the 
vaporization process (refs. 2 and 3) show an acoustic gain that is insufficient to over- 
come acoustic losses and cause instability. Acoustic oscillations with sinusoidal wave 
shapes were assumed in the previous studies; wave distortion was neglected. 

The importance of wave distortion in unstable combustion was first observed in 
studies with numerical models of the type developed by Priem and Guentert (ref. 4). 

These two-dimensional model studies of the combustor (both the combustion mechanism 
and the gas dynamics are modeled) show that a drop -vaporization process can support 
unstable combustion. An explanation of the result is not obvious, and the validity of the 
result could be questioned. However, examining the detailed numerical computations 
showed that wave distortion could be very important. Althou^ a sinusoidal disturbance 
is initially introduced, the pressure-time histories presented in references 4 and 5 show 
that this property is quickly lost and that the instability develops and grows into a dis- 
torted wave form. Acoustic gains extracted from the distorted -wave instability in refer- 
ence 5 are about three times larger than those obtained from the dynamic response of 
the same vaporization process to sinusoidal oscillations. These observations prompted 
the present study. 

A vaporization model where burning rate is simply related to some power of the drop 
Reynolds number is used. As such, the model generally characterizes convective heat- 
and mass -transfer processes, and the results are qualitatively applicable to a variety of 
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combustion mechanisms. This study, however, is specifically directed toward drop 
vaporization. 

The open -loop response of the vaporization process to a variety of wave shapes is 
analyzed. Approximate analytical expressions are derived for both the in-phase and out- 
of -phase response factors (components of the burning rate in-phase and out-of-phase 
with the pressure oscillation). Solutions for some specific conditions of wave distortion 
are presented, and the probable effect on unstable combustion is discussed. Exact so- 
lutions obtained by numerical techniques are also presented. They show the precision 
obtained by the theoretical approach. Some of these numerical results have previously 
appeared in reference 6. 


OPEN-LOOP MODEL 

The open -loop response of an assumed combustion process will be used to predict 
the acoustic amplification of the process within an unstable combustor. In such open- 
loop analyses the perturbation in burning rate caused by assumed acoustic oscillations 
are evaluated. The response of the process is obtained by comparing the output (burning 
rate) with the input (acoustic oscillations). 


Combustion Process 

Precise modeling of a specific combustion process is not intended or necessary for 
the purpose of this analysis. It will be shown that a combustion process where burning 
rate (mass or energy release rate per unit volume), simply given by 

w = C^(Re)“ = ( 1 ) 


is a functional relation particularly sensitive to wave distortion. (Symbols are defined in 
appendix A. ) 

A Reynolds number dependency (eq. (1)) characterizes the convective heat- and 
mass -transfer processes that dominate many combustion mechanisms. Among these 
are the erosive burning of solid propellants and drop vaporization in liquid rockets and jet 
engines. Drop vaporization is particularly dominated by such Reynolds number de- 
pendence, and it will be emphasized in this analysis. For drop vaporization the burning 
rate is usually proportional to the one-half power of the drop Reynolds number, m = 1/2. 

Priem and Guentert (ref. 4) used a vaporization model equivalent to that given by 
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equation (1) to obtain their useful correlations for nonlinear instability limits of liquid 
rockets. The approach they originally used will be followed. The gas viscosity p and 
the drop diameter D will be assumed to be constant. 

The relative velocity U in the Reynolds number (eq. (1)) introduces an important 
distortion-sensitive parameter. For an axially moving drop in a cylindrical chamber, 

U is the magnitude of the relative velocity vector given by 


Gas velocity oscillations in any of three vector directions affect the burning rate. 
Because of the greater interest in transverse -mode instability in rocket combustors, 
this analysis will be restricted to transverse -velocity oscillations. The axial -velocity 
difference |u, - u I is assumed to be constant. For simplicity, the transverse ve- 
locities are combined into a single transverse -velocity vector u^ giving 


U = AV 


1 + 



( 3 ) 


where 


AV = 



u 


2 

t 



+ u 


2 

r 


When a perturbed gas density is defined as p = P(1 + P’)> a burning rate normalized by 
the burning rate without oscillations w/wq can be expressed as 


where 


W = 


w 


w. 


= (1 + p’)"^ 




-,m/2 


( 4 ) 
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Equation (4) specifies the burning rate for any assumed transverse acoustic oscillation. 
As in reference 4, the process as modeled is insensitive to the frequency of the oscil- 
lation. 

Actually, the response of the vaporization process is frequency dependent (refs. 2 
and 3). Figure 1 shows a tjqjical response to sinusoidal oscillations. (Response is de- 
fined on p, 6. ) The peak value of the response denotes the most unstable condition. This 
peak value depends on the Reynolds number sensitivity of the process and is proportional 
to the exponential power of the drop Reynolds number. Oscillations in drop diameter 
reduce the response at low frequencies, and oscillations in drop temperature reduce it 
at high frequencies. Without these frequency -dependent properties, this analysis will 
only show how this frequency-response curve is modified by wave distortion. It will 
show how much the response, particularly the peak value, could be amplified or atten- 
uated by distortion waves. 


Acoustic Oscillations 

The acoustic oscillations assumed for this analysis are given by 


P' = ^ Pj^cos(no)t - cp^) 
n=l 


P' = Pj^cos(nwt - cp^ 
n=l 


u^ = ^ cUj^ cos(nojt - cp^- e^) (5) 

n=l 

The oscillating gas properties are expressed as Fourier series which can describe 
any wave shape. The phase angle cp^ specifies the phase relation between harmonic 
components. This harmonic phase angle is assumed to be identical for pressure and 
density, an approximation that is reasonably realistic. The gas velocity expression also 
contains this harmonic phase angle but is further modified by a velocity phase angle 0^^. 
This velocity phase angle specifies the velocity -pressure phase relations. When 0^^ is 
zero, the velocity and pressure are in-phase, and traveling wave properties are sim- 
ulated. For standing acoustic modes, the velocity and pressure are out of phase; a 
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9^ of 90° simulates this condition. Some acoustic modes exhibit neither pure traveling 
or standing properties. 


Response 


Imposing the acoustic oscillations expressed as harmonic series on the combustion 
process causes a perturbation in the burning rate that can also be expressed as a har- 
monic series. An important parameter obtained from a comparison of such burning- 
rate oscillations with the acoustic oscillations is the portion of the burning rate that is 
in phase with the pressure oscillations. Based on the Rayleigh criterion for heat driven 
waves, this in-phase component gages with the ability of the process to drive the acoustic 
oscillations. 

The in-phase component of the burning rate will be called the nonlinear in-phase 
response factor, It can be extracted from the burning rate and normalized by a 

correlation procedure (ref. 2) defined by 


(R 


nl 



( 6 ) 


For rocket combustors, in-phase response factors greater than some value between 0. 8 
and 1.0 will usually denote unstable combustion for sinusoidal oscillations; that is, the 
oscillations will grow with time. Whether the same critical values apply to the growth 
of distorted waves has not been rigorously established. Intuition and an examination of 
the distorted wave instability with the Priem-Guentert model indicate that the same crit- 
ical values apply. Without analytical proof at this time, a nonlinear in-phase response 
factor of 0. 8 to 1. 0 will be assumed to discriminate between stable and unstable 
combustion. 

The portion of the burning rate that is out-of -phase with pressure is also important 
in combustor stability studies with sinusoidal waves. This out-of -phase component 
largely effects only the frequency of an unstable combustor. However, the exact role 
and definition of an out-of -phase response factor for distorted waves also lacks analytical 
treatment. 

The nonlinear out-of-phase response factor defined by 
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( 7 ) 


r2n ^ 

j W’ / > Pj^sin(no>t - 

f = *U _ 

nl ~ 


/ 2ir 

(p 


O^darf: 


will be used. The definition is similar to that used for sinusoidal waves. This non- 
linear out-of -phase response factor together with the nonlinear in-phase response 
factor are considered to be an adequate characterization of the response of a pro- 
cess to distorted acoustic oscillations. 

Response factors can be used in the linear analysis of combustor stability to pre- 
cisely couple the combustion process to the gas dynamics. In linear analyses, where all 
oscillations have sinusoidal wave shapes, the relationship 

W* = ((R + iy)P’ (8) 


exactly specifies the coupling. 

When distorted-gas dynamic waves cause distorted buming-rate perturbations, a 
precise coupling relation has not been established. However, approximate relations 
can be formulated. This study is directed toward the approximation given by 

i>j^P j^Pl cos(a>t - (p^) + P 2 cos( 2 o)t - ^ 2 ) + • • • 

+ Pj^ cos(no)t - + residuals (9) 

With this approximation, the coupling is established by the two nonlinear response fac- 
tors and which will be discussed throughout this report. 

These nonlinear response factors are in fact weighted averages of harmonic re- 
sponse factors that can be specified for each harmonic component of the acoustic oscil- 
lations. For example, the nonlinear in-phase response factor given by equation (6) can 
be shown to be equivalent to an average of the harmonic response factors given by 


2 2 
Pj(Rl + P2<R.2 




^nl 


2 2 

Pi + P2 + 


( 10 ) 


+ p; 
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where the harmonic response factors are specified by 



Another approximate coupling relationship, formulated by conventional methods of 
harmonic ordering, weighs the harmonic response factors in a different manner. The 
relationship formulated by such methods is 

W* = ((R j + i>j)pj cos(o)t - <pj) + ((R 2 + cos(2o)t - + . . . 

+ ((Rj^ + cos(nwt - cp^ + residuals (12) 

With this approximation a series of harmonic response factors must be specified to 
establish the coupling between the gas dsmamics and the combustion process. (A dis- 
cussion of the harmonic response factors is given in appendix B.) 

There is no rigorous justification for the use of either of the coupling relationships 
presented, and others could also be formulated. The degree of precision obtained by 
any approximation will depend on the specific application and the method of handling the 
residuals. 

THEORETICAL SOLUTIONS 

The method used to obtain anal 3 dical solutions for the nonlinear response factors 
will be described. A general solution will be derived. The composition of the general 
solution and the reasons for the effect of distortion variables on response will be dis- 
cussed. 

Some additional simplifying assumptions will be made to describe the acoustic os- 
cillations. With these assumptions analytical solutions having greater utility are ob- 
tained. These solutions will be used in the RESULTS AND DISCUSSION section to des- 
cribe the effect of some specific variations in wave distortion on the response factor. 
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General Solution and Discussion 


Analytical solutions for the response factors were obtained by assuming exponential 
notations to simplify the mathematical procedures. The burning-rate expression 
(eq. (4)) combined with the assumed acoustic oscillations (eq. (5)) gives 


W = 


1 + 


OO 

Z 

n=l 


2 2 


m 


X 


1 + 


■ OO 

Z^-' 


.m/2 


i(nwt-<;Oj^-e^) _i(no)t-<p^-0^) 

2 


n=l 


(13) 


Two radical terms appear in equation (13). Each was separately expanded by a 
Taylor series. However, the squaring of the harmonic function in the radical term on 
the right was performed before the expansion. Both radical terms were expanded about 
their time invariant quantities, and only the first terms of the expansions were used. 

Multiplying the two expansions as indicated by equation (13) gives an expression for 
the burning rate of the form 


W = W(1 + W’) 


(15) 


The perturbation in burning rate W* can be arranged in harmonic series given by 


OO OO r .: 


i(nwt+o!_ J 


n, r' 


+ e 


-i(no;t+a!^ r)j 


(16) 


n=l r=l 


which is of the same form as the pressure perturbation 


OO 

-Z?[-' 

n=l 


(nwt-<pj^) -i(nwt-<p 


„)] 
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In equation (16) the A terms are generated by the cross product of the pressure 

n, r 

and velocity coefficients p and u . The a terms are sums and differences of 
the harmonic and velocity -pressure phase angles <p^ and 0^. 

The nonlinear in-phase response factor is obtained by performing the integra- 
tion indicated by equation (7), which involves the product W’P*. Performing the inte- 
gration gives 



(17) 


which can also be expressed as 


(R 


nl 


n=l r=l 



(18) 


A solution for where acoustic properties higher than the second order in har- 
monic content and cross products of harmonic coefficients higher than fourth order are 
neglected is given by 


' ^ 1^1 ^2^2 ^ 1 /j_Y 
. 2 ’ 2 2 ( AV I 

i Pj + ?2 ' ' 


PjUjU2Cos(c'>2 + ^2 " *^1^ 2 ~ + P2P2)^l‘^2 ^1 


— ' — ^ ^ 

cos 0^ + H 

^ a 


PlPfUj cos 20 J + — P2P2^2 ^^2 




(19) 


In this solution cp-^ = 0 has been assumed; that is, all phase angles have been refer- 
enced to the first harmonic pressure oscillations. 

Some qualitative interpretations of equation (19) are possible. The response is 
basically proportional to the exponent m on the Reynold’s number in the burning-rate 
expression. The solution consists of two major divisions: term 1 and terms 2 to 8. 
Term 1 is comparable to the solution obtained by linear (sinusoidal perturbation) analy- 
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sis. Its value is close to one and it is insensitive to variations in acoustic properties 
when P„ = ipn- 

The terms within the bracket are the primary contributions of wave distortion to 
the response. Their effect on the response factor increases with a decrease in the 
steady axial-velocity difference AV. All these terms depend on the phase relations used 
to describe wave distortion. Maximum values are obtained for zero phase angles (i. e. , 
velocity in-phase with pressure and no phase displacement between harmonic compo- 
nents). The response factor is amplified by wave distortion for such in-phase acoustic 
properties. Attentuation occurs for certain combinations of phase angles that give large 
negative values for the cosine functions. 

Terms 2 and 3 are the major contributions of wave distortion to the response be- 
cause of their lower order in cross products of the coefficients. Term 2 arises because 
the product of the first- and second -harmonic velocity components produces a first - 
harmonic component in the burning rate. This harmonic can be in-phase with the first - 
harmonic component of the pressure oscillation. In term 3 the square of the first har- 
monic terms of the velocity produces second -harmonic burning-rate components, which 
can correlate with the second-harmonic pressure component. Terms 4 to 6, although 
high order, can be interpreted in a similar manner. Although not so arranged, equa- 
tion (19) is composed of first- and second-harmonic responses and (fl 2 > which com- 
bine to give the nonlinear in-phase response The contributions to (R^^ and (R 2 can 

be recognized by the respective pressure coefficients Pj^ and P 2 appearing in terms 2 
to 6. (Appendix B gives solutions for (Rj and (R 2 -) 

Terms 7 and 8 can be considered to be normalizing factors. Term 7 arises from 
the numerator in equation (7) and is basically the square of the rms pressure oscillation. 
Term 8 is related to the mean value of the burning rate, where 



The mean burning rate and term 8 increase with an increase in the amplitude of the 
velocity oscillations, but is independent of the pressure amplitude. 


Specific Solution 

The general solution given by equation (19) can be expressed in a form that is more 
conveniently used to examine some specific oscillating conditions. Some interrelations 
between acoustic gas properties are assumed for this purpose. The gas density per- 
turbations are assumed to be proportional to the pressure perturbations such that 
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T^n = Pn' equation is usually valid only for small perturbations. Also, the co- 

efficients of the gas velocity expressions are assumed to be proportional to the pressure 
such that >u^ = This is also adopted from small perturbation theory. The pro- 

portionality constant is equal to one for weak traveling waves. 

With these restrictions and definitions and making specific application to drop va- 
porization where m = 1/2, equation (19) becomes 



Evaluations of equation (20) are used in the RESULT AND DBCUSSION to examine the 
effect of wave distortion on the response. 

A solution for the nonlinear out -of -phase response can be obtained by a method 

similar to that used for the (R„, solution. The result for rP^ = p , yu = J p , and 

m ‘ n '^n' ' n n^n 

m = 1/2, a solution comparable to equation (20), is 



This out-of-phase response is zero when the phase angles are zero, that is, when 
pressure -velocity and the harmonic components are all in-phase. At other phase angles 
it varies inversely with the relative axial gas velocity, AV. Most of the terms in equa- 
tion (21) are similar to those for in equation (20) except for a sine and cosine inter- 
change of the phase angle functions. This comparison of equations (20) and (21) shows 
that a change in the in-phase property of a burning-rate perturbation causes a com- 
plimentary change in the out-of -phase property. 


NUMERICAL SOLUTIONS 

The analytical solutions for the response factors derived in the previous section are 
approximate because (1) Taylor series expansions were used and (2) some higher order 
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cross products of coefficients were neglected. Exact solutions, which do not have these 
limitations, were obtained by numerical techniques for the conditions described in the 
RESULTS AND DISCUSSION section. These exact solutions, when compared with the 
analytical solutions, show the overall precision of the more readily evaluated analytical 
expressions. 

Numerical solutions were obtained by programming the burning rate, the acoustic 
oscillations, and the response expression for exact solutions on a digital computer. The 
solutions were restricted to the conditions that apply to the analjdical solutions given 
by equations (20) and (21). The expressions used are as follows: 

Burning rate - 


W = (1 + 


1 + 



ql/4 


Acoustic oscillations - 


P’ = Pl 


cos cut 


P2 

+ — cos(2cut - 

Pi 


(P 2 ) 



cos cut + 


P2 

— cos (2 cut - 

Pi 


(P 2 ) 


Ut = 


^IPl 

r 


Jo P5 

cos(cut - 6^ + — — cos(2cut - (po “ 

h Pi 


Numerical integration of equations (7) and (10) gave the response factors and 


RESULTS AND DISCUSSION 

The results to be shown and discussed will be those for a combustion process where 
w « (Re)^'^^. This process applies most directly to drop vaporization in rocket com- 
bustors. The variables which describe the distortion of acoustic oscillation will be ex- 
amined. Generally, the independent effect of wave distortion variables on the response 
of the process will be shown. The purpose of the discussion is to demonstrate the im- 
portance of wave distortion in unstable combustion. 
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Harmonic Content 


The effect of varying the amount of second harmonic content on the nonlinear in- 
phase response is shown in figure 2. The effect on pressure wave shape is also 
shown. For the conditions of figure 2 all the phase angles are zero; that is, there is no 
phase displacement between harmonic components and a velocity in-phase with pressure 
and density. The relative axial Mach number term AV/cJj is 0.02, which is a rela- 
tive drop velocity in rocket combustors of about 30 meters per second (100 ft/sec). 

The in-phase response exhibits optimum properties as a function of both pressure 
amplitude and second harmonic contents. Maximum response occurs at an amplitude of 
about 0. 02 for a second -harmonic to first -harmonic amplitude ratio of about 0. 8. One 
of the reasons for the low response at higher amplitudes is the increase in mean burning 
rate with amplitude. Increasing the mean decreases the relative size of a perturbation 
and reduces the response (see appendix C for mean values) . 

The significance of the response values shown in figure 2 can be gaged by a compari- 
son with the P2^^1 ~ curve in figure 2. The P 2 /P 1 = 0 curve is the response to sinus- 
oidal waves and is shown to be relatively insensitive to pressure amplitude. The re- 
sponse is about l/2y (the value obtained by linear analysis), which is too low to cause in- 
stability in rocket combustors. Thus, distorting the acoustic oscillations with second- 
harmonic content can amplify the response by a factor of 10 above this linear value. The 
response is much larger than that needed for instability (i. e. , 61^^^ is greater than 0. 8 
to 1.0.) The contribution of the harmonic response factors and (R 2 to the overall 
nonlinear response is discussed in appendix B. 

For the conditions shown in figure 2, maximum amplification occurs at the rela- 
tively low oscillation pressure amplitude of 0.02. This amplitude is comparable to the 
usual combustion noise in rocket engines. It is possible that combustion instability could 
develop from distorted acoustic oscillations buried within the usual noise spectrum of 
the engine. 

A comparison of the analytical and numerical results in figure 2 shows that the analy- 
tical solutions overpredict the response by about 10 to 20 percent. The functional form 
of the solution, however, appears to be correct. 

Some types of wave distortion are not adequately characterized by only two-harmonic 
components as discussed thus far in this report. This is particularly true for high- 
amplitude combustion oscillations where the wave shape becomes highly distorted. A 
high -amplitude wave shape of special interest in rocket -engine instability is that for the 
traveling transverse acoustic mode. High amplitudes for this mode are frequentily ob- 
served during instability with large rocket engines. 

The pressure wave for strong traveling acoustic modes has a characteristic shape 
which was experimentally examined and described in reference 7. The harmonic content 
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of increasingly higher orders was needed to characterize wave shapes for increasing 
amplitudes. The following relation approximated the behavior 




cos nwt 


( 22 ) 


where the maximum to minimum and mean pressure amplitudes are given by 


p» _ p» . = 

max mm 


2Pi 


1 - Pi 


(23) 


and 


Pi 


rms 


V2(l - Pi> 


(24) 


The following table shows the harmonic content and amplitudes: 


Pi 

P 2 /P 1 

P3/P1 

P4/P1 

P' - P’ . 

max min 

P' 

rms 

0.01 

0.01 

0.0001 

0.000001 

0.02 

0.0071 

. 1 

. 1 

.01 

.001 

.201 

.071 

. 3 

. 3 

.09 

.027 

.66 

. 165 

. 5 

• 5 

.25 

. 125 

1. 33 

.41 


Numerical solutions for the response factors were obtained for these waves. As 
previously, the assumptions that yP^ = p^^, yu^ = and = 1.0 were made. The 

results and typical pressure wave shapes are shown in figure 3. The response increases 
with amplitude and assymptotically reaches a maximum value of about 1. 08. (The var- 
iation with relative velocity will be discussed later. ) The response does not reach the 
peak value shown in figure 2 because the harmonic content for these waves is very small 
at the pressure amplitude of 0. 02 where the peaks occur in figure 2. At this amplitude, 
P 2 /P 1 is only 0. 02 for these traveling mode waves. Also, the response does not fall at 
high amplitudes as in figure 2 because higher order harmonic components contribute to 
the response in this region (see appendix B). The maximum response of 1.08 (about 
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2. 5 times the linear value) could cause unstable combustion. The harmonic distortion 
which grows with amplitude, therefore, is a potential cause of instability in rocket en- 
gines. 

The results in figure 3 provide some insight into the causes for trigerable insta- 
bility observed both experimentally and with the numerical model studies of Priem and 
Guentert (ref. 4). When disturbances of increasing amplitude are introduced in a com- 
bustor the harmonic content can be expected to grow with amplitude. At some amplitude 
there is sufficient harmonic content to give an in-phase response that can cause insta- 
bility greater than 0. 8 to 1. 0). It appears that Priem and Guentert’s numerical 
model studies depended on this amplification caused by wave distortion. It is probable 
that instability in rocket engines can be triggered by bomb pulses when they initiate os- 
cillations with sufficient harmonic content. 

Figure 3 also shows analytical results obtained from equation (20). The substitution 

2 2 
P 2 = (pj) (from eq. (22)) was used and terms of a higher order than pj were neglected. 

This, together with the assumptions already stated, gives the result 



This solution good to second order in harmonic content and amplitude coefficient 
agrees favorably with the exact solution for multiharmonic wave shapes. Some addi- 
tional evaluations showed that deductions with regard to phase relations can also be 
made for these multiharmonic wave shapes on retaining angles in the derivation of 
equation (25). Also, a solution for / for strong traveling waves can be deduced 
from equation (21). 

Pressure-Velocity and Harmonic Phase Relations 

The amplification caused by wave distortion depends on both the phase relation be- 
tween the first- and second-harmonic components cp^ and the phase relations between 
pressure and velocity 0^. The solutions for the response factors show that phase angle 
effects are not separable in either equation (20) for or equation (21) for Many 

interacting effects are possible. The discussion of the phase relations will be limited to 
some of their gross effects. 

The independent effect of the harmonic phase angle on the in-phase response factor 
(Rj^^ is shown in figures 4(a) and (b). The conditions are comparable to figure 2, but 
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only a second to first harmonic amplitude ratio of 0. 2 is considered. The pressure- 
velocity phase angle is zero (0^ = 02 = 0°)* The in-phase response is shown to decrease 
for any phase shift between harmonic components. When the phase shift is 90°, the re- 
sponse is reduced to near the linear value of l/2y. Phase shifts of more than 90° can 
give negative values of the in-phase response. Negative responses imply highly stable 
combustion. 

Application of these harmonic phase-shift results to combustion instability is quali- 
tative because the phase -shift properties of actual waves have received little attention. 
Steep -front or shock-like waves contain a second harmonic component which approaches 
a phase shift of 90°. Such waves would not be as easily driven by the combustion pro- 
cess as are the strong symmetrical wave shapes shown in figure 3 where there is no 
harmonic phase shift. Phase shifts beyond 90° seem improbable except for some types 
of multimode resonance. Any method of producing large harmonic phase shifts, how- 
ever, should contribute to combustion stability. 

The effect of harmonic phase shift on the out-of-phase response is shown in fig- 
ures 4(c) and (d). Any phase shift between harmonic components causes an out-of -phase 
response. When the phase shift is positive, the response is positive; when the 

phase-shift is negative, J is negative. The result is unusual because out-of-phase 
responses in linear analyses appear only for time -dependent processes. The signifi- 
cance of these out -of -phase responses from wave distortion remains to be established. 

The independent effect of the velocity -pressure phase relation on the in-phase re- 
sponse is shown in figures 5(a) and (b). Any phase shift between pressure and velocity 
reduces the response, but the effect is much smaller than for a harmonic phase shift. 
The maximum reduction occurs for a 90° phase shift, but the response remains above 
the linear value of l/2y. A velocity out-of-phase with pressure (0 = 90°) characterizes 
standing mode resonance. The result implies that the standing mode resonance is not 
as easily driven as is the traveling mode resonance which exhibits in-phase velocity- 
pressure properties. 

Any phase shift between pressure and velocity also causes an out-of-phase response 
*^nZ ^° This is shown in figures 5(c) and (d). The response is zero at 90° 

increments rather than the 180° increments found for the velocity-pressure angles in 
figure 4. 

One condition of practical interest is an increase in the harmonic phase angle when 
the velocity -pressure phase angle is 90°. This simulates the growth of shock-like 
waves in standing modes. The velocity is assumed to be out -of -phase with pressure 
(0j^ = 02 = 90°) because of the standing mode property. An increase in the phase angle 
between the harmonic components provides the wave steepening property. Figure 6 
shows the effect on both and The effect of the harmonic phase angle on 

is similar to that when the pressure -velocity phase angle is zero (fig. 5). A harmonic 
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phase shift of 90° (a limiting condition for shock-like waves) reduces the response to its 
linear value of l/2y. Higher harmonic components should also be introduced for a 
complete solution for shock -like waves. It is concluded, however, that steepening the 
wave reduces the response, which in turn acts to limit the pressure amplitude; that is, 
a wave steepens when it grows in amplitude until the acoustic gains represented by the 
in-phase response match or equal the acoustic losses. 


Relative Axial Velocity 

The velocity sensitivity of a Reynolds number dependent process is the primary 
cause of the large changes in response with wave distortion. For drop vaporization this 
velocity sensitivity is controlled by the steady relative axial -drop velocity or Mach num- 
ber AV/c. When this relative velocity is large, the vaporization rate is large. This 
makes the process relatively insensitive to velocity oscillations and any distortion of 
such oscillations. 

The effect of the steady relative axial velocity on the in-phase response is shown in 
figure 7. For very large relative velocities, the response approaches the value of l/2y 
obtained by linear analysis. For the conditions used in figure 7, the amplification caused 
by wave distortion becomes appreciable only for relative axial Mach numbers less than 
about 0. 1 (about 150 m/sec or 500 ft/sec in rocket combustors). However, Mach num- 
bers below 0. 1 typify drop vaporization in many rocket combustors, and a significant 
increase in response due to distortion can be expected. 

The effect of the steady axial velocity is also shown in figure 3 which shows the in- 
phase response for strong traveling waves with multiple harmonic components. For 
these waves a decrease in the relative velocity decreases the pressure amplitude needed 
to obtain large response factors. 

The relative axial velocity was one variable examined by Priem and Guentert (ref. 4) 
in their numerical model studies of instability. They showed that the size of the dis- 
turbance needed to start instability decreased with a decrease in relative velocity. The 
same result is implied by figure 3. The pressure amplitudes needed for instability 
^'^nl 0. 8 to 1. 0) decrease with a decrease in relative velocity. 

The trend for large-scale liquid rockets is toward smaller contraction ratio designs 
which lead to higher relative axial Mach numbers. For example, an average AV/C of 
about 0. 15 characterizes drop vaporization in a relatively small contraction ratio com- 
bustor of 1.8. These smaller contraction ratio combustors should be more stable, but 
stability is not assured. As shown in figure 3, an increase in relative Mach number 
merely increases the pressure disturbance that the process can tolerate and still remain 
stable. 
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CONCLUDING REMARKS 


This analysis using a simplified vaporization process has shown wave distortion to 
be an important variable in combustion stability. It can amplify the response of the pro- 
cess beyond that needed for unstable combustion and significantly affect stability pre- 
dictions. This effect of wave distortion on stability should persist for a more precise 
vaporization model or for any combustion mechanism where Reynolds number is a dom- 
inant variable affecting combustion rate. It appears necessary to direct more attention 
to the causes and effects of wave distortion to obtain a complete understanding of the un- 
stable combustion system. 

Various deductions about combustor behavior can be made from this study and com - 
pared with experimental evidence. But such comparisons must be qualitative for sev- 
eral reasons. One, the analysis is incomplete. Only a simplified combustion process 
was examined, and no overall systems stability analysis was made. Two, conflicting 
experimental evidence is often found, or the independent effect of a particular param- 
eter is not conclusively established. Recognizing these limitations, the following are 
of interest: 

(JO Traveling mode resonance is m ore freq uen tly enc oun tered than standing mode 
resonance . - Acoustic gains are larger for the in-phase velocity-pressure properties of 
traveling waves than the out-of -phase properties of standing waves. 

(2) Pressure ampliti^es for traveling mo^de instability are us ually larger than those 
for longitudinal standing mo de instability . - Acoustic gains are reduced by the steepening 
of the wave shapes, which develops naturally in longitudinal modes at relatively low am- 
plitudes. 

(3) Disturbances above some finite amplitude can trigger unst able combustion . - 
Acoustic gains increase with the increase in harmonic distortion, which usually grows 
with amplitude. 

(4) Large contraction ratio combustors are more prone to transverse mode instability 
than with low contraction ratios . - Acoustic gains from distorted waves are larger for 
the small relative axial drop velocities usually found in large contraction ratio combus- 
tors. 

(5) Evidence of multimode resonance is common during the initial growth period of 
spontaneous instability . - Acoustic gains are large at low amplitudes for waves distorted 
by two harmonic components. 

(6) Acoustic liners with broad-band tuning are most effective in suppressing insta- 
bility . - A reduction in harmonic content (from broad-band tuning) may be more effective 
in reducing acoustic gains than reducing only the fundamental oscillation. 

(7) The symmetrical wave shapes of traveling transverse mode instability tend to 
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steepen as they appr oach a limiting amplitude. - Acoustic gains decrease with wave 
steepening and eventually act to limit the amplitude. 


SUMMARY OF RESULTS 

The response of a combustion process to acoustic oscillations distorted by harmonic 
content was analyzed. A combustion process where burning rate is proportional to an 
exponential power of Reynolds number was assumed. Specific application was made to 
drop vaporization in rocket combustors where the exponential power is 1/2. The follow- 
ing results were obtained: 

1. Anal 5 dical solutions were obtained for the in-phase and out-of -phase response 
factors (components of the burning rate in-phase and out -of -phase with the pressure os- 
cillation, respectively). The solutions show the effects on the response factors of the 
harmonic amplitudes and phase relations, of the acoustic velocity amplitudes and 
velocity -pressure phase relation, and of a steady relative axial velocity. The analytical 
solutions show good agreement with exact numerical solutions. 

2. Harmonic content can substantially increase the in-phase response factor. The 
addition of a second harmonic component to the acoustic oscillation can increase the re- 
sponse by an order of magnitude above that for sinusoidal oscillations. Higher order 
harmonic distortion (the type that normally grows with amplitude) can increase the in- 
phase response to 2. 5 times the sinusoidal value. 

3. The phase relation between harmonic components and the phase relation between 
velocity and pressure affect both the in-phase and out-of -phase response factors. The 
in-phase response factor is a maximum when both the harmonic components and the 
velocity -pressure are in-phase. No out-of -phase response exists for this condition. The 
in-phase response factor is less than the maximum when either the harmonic components 
or the velocity-pressure are phase shifted. With phase shifts, the out-of -phase response 
factor usually has a finite value. The effects of the harmonic component phase relation 
are larger than those for the velocity -pressure phase relation. 

4. The relative axial velocity (of a burning drop in a combustor) controls the sensi- 
tivity of the process to wave distortion. Distortion effects on both response factors are 
significant when the relative axial Mach number is less than about 0. 1. 

5. Wave distortion is a crucial property affecting the response of the drop- 
vaporization process to acoustic oscillation. Changes in distortion can amplify or sup- 
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press the in-phase response factor and substantially affect the stability prediction of 
rocket combustors. Wave distortion should be carefully considered in the design of 
stable combustors. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 6, 1971, 

128-31. 
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APPENDIX A 


SYMBOLS 


Vr 

harmonic coefficient for burn- 
ing rate perturbations 

Uz,Ug,Ur 

gas velocity (axial, tangential, 
and radial) 

Cl 

proportionality constant 

"t 

2 

transverse velocity, (u^ + 
2 1/2 ^ 

c 

speed of sound 


D 

drop diameter 


r' 

relative velocity of drop. 


out -of -phase response factor 


I"Z - "z 1 

(harmonic and nonlinear) 

W 

dimensionless burning rate 

Jn 

proportionality constant in 
>^n = JnPn 

w 

burning rate (energy release 
per unit time and volume) 

m 

Reynolds number exponent 

"^0 

initial burning rate without 

n 

harmonic order 


oscillations 

P 

pressure 

O' 

n, r 

harmonic phase angle for burn- 

^rms 

rms pressure amplitude, 


ing rate perturbations 



y 

ratio of specific heats 



^n 

harmonic phase angle for 
velocity -pressure relation 

Pn 

harmonic coefficient for pres- 

P 

gas density 

sure perturbation 

^n 

harmonic coefficients for 

Re 

Reynolds number 

density perturbation 

^n’^nZ 

in-phase response factor 
(harmonic and nonlinear) 


harmonic phase angle for 
pressure perturbations 

t 

time 

M 

gas viscosity 

U 

magnitude of relative drop 

CO 

frequency 


velocity vector 

f 

dimensionless perturbations. 


axial velocity of drop 


II 


harmonic coefficient for trans- 
verse Mach number pertur- 
bations 


mean values 
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APPENDIX B 


HARMONIC RESPONSE FACTORS 


In this appendix harmonic response factors are presented for some representative 
conditions to show the magnitude of their contributions to the nonlinear response factor . 

The expressions developed for the relation between the nonlinear and harmonic re- 
sponse factors (eqs. (8) and (11)) are 


OO 


(R 7 = 5-Jl_ 

VLL OO 


2 (Pn)' 


n=l 


OO 


^nl ■ oo 


Z)<Pn)‘ 


n=l 


(Bl) 


(B2) 


The nonlinear response factors are shown to be weighted averages of the harmonic re- 
sponse factors. 

Both analytical and numerical solutions for the harmonic response factors were ob- 
tained. Analytical expressions for the first- and second-harmonic response factors 
were obtained by a method similar to that used for the nonlinear response factors. The 
expressions for yp^ = p^^, Jj = J 2 = 1.0, and m = 1/2 (comparable to eqs. (20) and (21) 
for the nonlinear response factors) are 
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Figures 8(a) and (b) show the in-phase response factors (R^, (R 2 > 3-nd for the 
highest value curve in figure 2. For this condition the harmonic components and the 
velocity -pressure are in-phase, (AV/cJ^) = 0. 02 and V 2 /V 1 = 0,8. Both the analytical 
and numerical solutions show (R^ to be larger than and (R 2 smaller. For this 

case of a relatively large amount of second -harmonic content, the weighted averaging 
gives a nonlinear response factor that is smaller than the first harmonic response factor. 

The relation between response factors are reversed when the second harmonic con- 
tent is reduced. This is shown in figures 8(c) and (d) for P 2 /Pj^ = 0-2. The second- 
harmonic response factor is very large for this condition. This makes the nonlinear 
response larger than the first-harmonic response. The second-harmonic response is 
large because the process considered is a strong second -harmonic generator or fre- 
quency doubler (see appendix C). When the second -harmonic content of the acoustic os- 
cillation is small, there is a strong correlation between the burning rate and this second- 
harmonic content. Increasing the second -harmonic content of the acoustic oscillation 
does not cause a proportional increase in the second -harmonic content of the burning 
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rate. This reduces the correlation and the second-harmonic response factor as ob- 
served in figures 8(c) and (d). 

The out -of -phase response factors exhibit a behavior similar to that for the in-phase 
response factors. An example of one condition is shown in figure 9. This is comparable 
to figures 8(c) and (d) in that the second-harmonic response is much larger than the 
first-harmonic response. The weighted averaging gives a nonlinear response that is 
slightly higher than the first harmonic response. 

Higher order response factors contribute to the nonlinear response for waves with 
multiple harmonic content. This is shown by the numerical solutions in figure 10 for 
strong traveling waves. As the pressure amplitude increases, higher order responses 
contribute to the nonlinear response and a constant value of nonlinear response is main- 
tained at high amplitudes. 


APPENDIX C 


INSTANTANEOUS AND MEAN BURNING RATES 

The burning rate, pressure, and acoustic velocity oscillations for a typical condition 
considered in this study are shown in figure 11. The burning-rate curve shows a large 
second -harmonic content even though the second-harmonic content of the acoustic oscil- 
lation is relatively small (P 2 /P 1 = 0.2). The process is basically a harmonic generator 
or frequency doubler. This causes the burning rate to be highly distorted for any type of 
disturbance. When the pressure and velocity are not in phase (as in fig. 11), the burning 
rate curve on the average is not in phase with the pressure oscillation. Such burning 
rates have both in-phase and out-of-phase components with respect to the pressure oscil- 
lation and the respective response factors and assume finite values. 

The mean burning rate of a drop vaporization or any Reynolds number dependent pro- 
cess generally increases in the presence of acoustic oscillations. 

The analytical expression for the mean burning rate for drop vaporization is given by 



For yu^ = Jj^Pj^ and yp^ = p^^, equation (Cl) becomes 



The mean burning rate is primarily dependent on the relative axial velocity and pres- 
sure amplitude although large second-harmonic contents also affect the mean value. 

The effect of relative velocity and amplitude is shown in figure 12 for P 2 /Pj^ = 0.2, 

Jf = J 2 = 1. 0, and y = 1. 2. The mean burning rate can more than double with an in- 
crease in pressure amplitude, and at low relative velocities the mean can increase by a 
factor of 5 or more. These increases in the mean with amplitude reduce the relative 
size of a given oscillation in burning rate and in part contribute to the decrease in re- 
sponse with amplitude shown in figures 2 and 4 to 7. 
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In-phase response factor, 



Characteristic frequency, dimensionless 

Figure 1. - Typical response of heptane drop vaporization to sinusoidal os- 
cillations (ref. 1). 


28 



Nonlinear in-phase response, 


Harmonic 
amplitude ratio, 
P2/P1 


0 



(a) Pressure wave shapes; pressure amplitude, constant. 



(c) Numerical solution. 


Figure 2. - Effect of amount of harmonic content on response and 
wave shape. Harmonic phase angle, (p 2 = 0°; velocity-pressure 
phase angles, = 0o = 0°; relative axial Mach number, 

AVfcJ] == 0.02; specific heat ratio, 1.2; proportionality constant, 
J2/Ji= 1.0. 


Nonlinear in-phase response, 


First order 
harmonic 

for pressure Pressure 
perturbation, amplitude, 

Pi Prms 



(a) Pressure wave shapes. 



Pressure amplitude, 


L i II II I I I I 

0 .01 .02 .05 .1 .2 .5125 

Maximum to minimum pressure amplitude, P[^gj^ - PJ^jp 

(c) Numerical solution. 

Figure 3. - Response and wave shapes for waves in strong traveling 
transverse acoustic modes. 
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Nonlinear out-of-phase response, Nonlinear in phase response, 


Harmonic 
phase angle, 


^ 2 ' 

deg 




(c) Analytical solution for 



(dl Numerical solution for 

Figure 4. - Effect of harmonic phase angle on response. Relative 
axial Mach number, AV/cJj = 0.02; specific heat ratio, y = 1.2; 
proportionality constant, harmonic ampiitude ratio, 

P2/P1 = 0- 2. 


Out-of-phase response, In-phase response, 


Velocity -pressure 
phase angle, 

8 , 



Figure 5. - Effect of velocity-pressure phase angle on response. 
Harmonic phase angle, (p 2 = 0°; relative axial Mach number, 
AV/cJ] = 0.02; specific heat ratio, y = 1. 2; proportionality 
constant, J 2 /Ji= 1 . 0 ; harmonic amplitude ratio, P 2 /Pi= 0 . 2. 


Out-of -phase response, In-phase response, 


Harmonic 
phase angle, 


deg 




Figure 6. - Effect of harmonic phase angle on response for an out-of- 
phase velocity-pressure. Velocity-pressure phase angles, = Go = 
90°; relative axial Mach number, AV/cJj= 0.02; specific heat ratio, 
y= 1.2; proportionality constant, harmonic amplitude 

ratio, P 2 /P 1 0. 2. 





response 

factors 



(a) Analytical solution for harmonic amplitude ratio P 2 /P 1 = 0.8. 


0 



(b) Numerical solution for harmonic amplitude ratio P 2 /P 1 = 0.8. 



(c) Analytical solution for harmonic amplitude ratio P 2 /P 1 = 0. 2. 



Pressure amplitude, Pj-ms 


(d) Numerical solution for harmonic amplitude ratio P 2 /px = 0.2. 

Figure 8. - Comparison of harmonic and nonlinear in-phase responses. 
Harmonic phase angle, ^ 0 “ 0°; velocity -pressure phase angles, 

Oj = 02= 0°; relative axialMach number, AV/cJi = 0.02; specific 
heat ratio, y= 1. 2; proportionality constant, J 2 /Ji= 1*0. 






f// Increments contributed by / 
y/. harmonic responses y 


Pressure amplitude, ^ms 

Figure 10. - Contribution of harmonic responses to nonlinear response for strong traveling waves. Numerical 
solution. 


.06r— 1.4, 




y%\ \p 


Figure 11. - Numerical evaluations of burning rate, pressure, and acoustic particle velocity 
perturbations. Harmonic phase angle, velocity-pressure phase angles, 0i = 02 

= 45°; relative axial Mach number, AV/cJj = 0.02; specific heatVatio, y = 1.2; proportionality 
constant, ^2lh^ harmonic amplitude ratio, P 2 /P 1 = 0.2; harmonic coefficient, p^= 0.3. 



Figure 12. - Mean burning rates. Harmonic phase angle, (p 2 "" \ velocity-pressure phase angles; 0^ = 02=O°; 

specific heat rates, y = 1. 2; proportionality constant, = 1.0; harmonic amplitude ratio, p^/pj = 0. 2. 
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